DUAL SCATTERING CHANNEL SCHEMES AND 
TRANSMISSION LINE METHODS 

STEFFEN HEIN 

Abstract. Dual scattering channel (DSC) schemes naturally extend 

the transmission line matrix (TLM) numerical method beyond the lines 

^|H ' set up by P.B. Johns and coworkers. Conceptually, DSC schemes re- 

("^ , tain from TLM the typical formal splitting of the computed fields into 

f^ ' incident and reflected components, interchanged and scattered between 

Cn ' adjacent cells, without yet a need for transmission lines. They gener- 

(-H . ally admit nontrivial cell interface scattering and so release from a set 

r-{ ' of modeling limitations to that the classical TLM method is subject. 

The elements characterizing DSC schemes are outlined and illustrated 
(^ ■ by coupling within a non-orthogonal mesh a model of heat propagation 

fT^ ' to a lossy Maxwell field. 

Westerham on February f , 2008 

Life is not an illogicality; yet it is a trap for logicians. 

It looks just a little more mathematical and regular than it is. 

Gilbert Keith Chesterton 

1. Introduction 
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Chesterton's trap lurs where theory encouters real life - or other theories. 
Mathematicians apply their methods, and sometimes export them to fields 
VJQ I where they are inadequate; just as physicists or microwave engineers intro- 

duce proper concepts into mathematics, with sometimes creative effects or 
walking into a trap - two issues that are not absent in the history of the 
(^ I transmission line matrix (TLM) numerical method. 

^" ■ In the early 1970s Peter B. Johns f , prominent professor of electrical engi- 

neering at the University of Nottingham, had with his coworkers the brilliant 
1-^ ' idea to apply transmission line network modeling techniques to the numeri- 

cs . cal solution of Maxwell's equations |.IoBe| iJoT] . There are in fact appealing 

H ! features of the TLM method, such as unconditional stability |,Io2j . high flex- 

ibility due to novel stub-loading (deflection) techniques |Ch|IHoe|lHelj . and 
easy implementation, which made the method attractive to users. TLM 
K> I approximation to wave propagation, diffusion, and transport phenomena 

j^ ■ became so a subject of assiduous study and publication [Tlml-3]. 

Familiarity with the transmission line picture and the well known scat- 
tering concept certainly fostered the acceptance of the TLM method in the 
microwave engineering community. However, the primary interest so turned 
of course on applications in their own discipline, rather than onto the inner 
structure of the TLM algorithm as an object of mathematical analysis. The 
underlying framework remained firmly linked to 'real' transmission lines, 
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indeed too firmly to allow for substantial evolution of the TLM method be- 
yond the original field of application. The algorithm could not unfold its 
entire potentiality. 

The history of the TLM method is best displayed in the proceedings of 
three conferences that exclusively focussed on the subject [Tlml-3]. Also, 
the monographs of Christopoulos |Ch| and de Cogan |dCj reproduce the 
original ideas and introduce into classical applications. Rebel |Rej gives a 
fairly complete survey over the state of the art by the year 2000 - which is 
more or less still the actual state. 

Dual scattering channel (DSC) schemes result from an incisive revision 
of the TLM paradigm. They are generalized transmission line methods in 
arising from a twofold abstraction: Firstly, the scattering channel concept 
underlying TLM is redefined in terms of paired distributions. (Character- 
istic impedances are thus neither needed, nor in general defined, e.g.) In 
the second place, nontrivial cell interface scattering is admitted during the 
connection step of iteration thus taking advantage of the intrinsic duality in 
the connection-refiection cycle of the algorithm. 

The new framework bypasses a set of modeling limitations imposed by 
transmission lines while it preserves the major advantages of the TLM 
method. In particular, the convolution type updating scheme and Johns' 
two step connection-reflection cycle are essentially retained. DSC schemes 
remain unconditionally stable under quite general circumstances, and they 
allow for process parallelization and deflection - the latter of which general- 



izes the stub loading feature of the classical TLM method [Hell IHe2j . 

DSC schemes are conceptually simple and natural. They are a coherent 
response to space time discretization in that meshing virtually imposes an 
artificial cell-boundary duality. A detailed technical introduction into DSC 
schemes has recently been given in jHe3j . The present paper offers a con- 
densed desription of the framework with emphasis on the typical contraction 
properties that ensure stability. 

Spinner's implementation of a DSC heat propagation scheme in non-or- 
thogonal mesh, wherein heat sources are coupled to a lossy Maxwell field, 
illustrates the approach. 

2. Causality 

Let £/ be a real or complex linear space and / a totally ordered set ( e.g. 
/ € {N , Z , M} ; intervals are then naturally defined in / by means of the 
order relation. Of course, we shall think with / of a discrete or continuous 
time domain.) Also, let £ C XL^ be a set of functions such that / G £ implies 
Xs<ti^)fi^) ^ ^' f°i' every t £ I , where Xp{s) denotes the characteristic 
function of property P (which is 1 if s shares that property and else). 
Sometimes we write simply [/]<( in the place of Xs<ti^) fi^)- 

Definition 1. A function F : £, ^ 8, is called causal, iff for every / £ £ 

""""^^^^ Ff{t) = F[x^^J{s)]{t) . 

Such functions are also called (causal) propagators. 
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Remark. In some respect, causal functions generalize lower triangular ma- 
trices or integral operators such as 

ft 



Ff{t) = f K{t- s)f{s)ds 



for instance with a Green's function kernel K. 

Note that in general F : f >-^ Ff needs not to be linear. 

Typically / represents a state evolving in time (i.e. a process). Then 
causality of F spells that F f [t) depends on the history of / only up to 
present time t . 

The proof of the following is easy and left to the reader. 

Proposition 1. 

For every t (z I ; f,g G £ and causal functions F, G : £. ^> S, 
(i) Ff{s) = F [ x,^t (r) f {r)]{s) , for every s < t . 

(ii) // f (s) = g (s) for every s ^ t , 

then Ff (s) = Fg (s) for every s ^ t . 

(iii) The product of causal functions 
FG : t — ^ £ 

/ ^ FGf := F[G[f]] 

is again causal. In fact, if £. is a linear space, then the causal func- 
tions over £ form an algebra. 

3. DSC Processes 

Just as the TLM algorithm, DSC schemes operate on a linear space T of 
so-called propagating fields , which is a product of (real or complex) normed 
spaces 

cf. |He3j . The two factors are (loosely) named the incident and outgoing 
subspace of J". They are isomorphic in that there exists a canonical involu- 
tary isomorphism of normed spaces 

nb : T ^ 7 

which is commonly called the node-boundary map. Hence, there exists a 
space [L , ||...|| ) such that 

rp ^^ IT) C^ ( r II II "l 

^ in — -^ out — V"^' lr"ll/ 

in the sense of isomorpy of normed linear spaces. In the same meaning 

'^ = (i:M|...|r) , 

e.g. with norm ||(a,6)||~ := \/||«|P + ||&|P ; a, 6 S XL 

(or any equivalent norm). 

DSC and TLM algorithms then follow a two-step iteration cycle in work- 
ing with alternate application of a connection and reflection map 

C : J X J Q^( > J in Ji : J X y^^ > J out 

{t, Zoutit - Kt)'^^q) ^ Zinit) , {t, Zin{t - Kt)'^^q) ^ Zout{t) 
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which respectively update the propagating fields at even and odd inte- 
ger multiples of half the time step, i.e. onI:= {kr \ k £ Z} and J : = 
{ {2k + 1)t/2 I A; G Z }. To these maps the following functions F^ , F^ are 
associated in a one-to-one correspondence 



c 



f ^ F^f with F^fit) := nboe{t,f: 



and 



F^ : i:-^ ^ £ 



J 

9 ^ Fr9 with F^g{t) := nboO{{t, g] 



Proposition 2. 

(i) F^ and F^ are causal on L^ and L"^ , respectively. 
(ii) For every r,s £ J and Tq : f[t) i-^ /(t -|- q) the shift operator on 

if r + s ^ , then < J, ^ „** is causal on < „ r 

The first statement is validated by observing that the reflection and con- 
nection maps !Jl and C are functions of back in time running sequences of 
propagating fields, cf. Q, and are themselves causal in this sense |He3j . The 
second statement is then a directly verified. □ 

For later use we fix straightaway: 

Definition 2. For any incident function e G iL^ x { } supported on a 
finite interval [ , Nt ) C / , the DSC process excited with e and generated 
by 01 and C is the unique function z = {zm , Zout ) G ( ■iJ^ y^"^ such that 
z (t) = for t ^ and recursively for < t G / U J 

(2) ....,■ ^ ,{^^nit),T_^^F^T_r[e + Zin]{t)) if t G / 

Thus, in particular, Zin and Zout 'switch' respectively at even and odd 
integer multiples of half the time step r. 

4. Passivity 

Algorithm stability prevents the computational process from piling up 
to infinity (of course, it does not yet imply convergence, or consistence of 
the algorithm). TLM models are unconditionally stable in that they are 
equivalent to passive linear transmission line networks |Jo2j . DSC schemes, 
in not using lines, need a more general characterization which here is given 
in terms of a-passive causal functions. 

Under the same assumtions for iL and / as in section |21 let ||...|| be a 
norm on L and a G K a continuous non-negative real functional. 

Definition 3. 

(i) We call a process g : I ^> H stable, iff it is uniformly bounded on / 

(i.e. there exists b G M+ such that |1 5 (t) || < b for every t £ I ). 
(ii) The functional a : £ — > M is named a (de) limiting functional, iff 
there exist any non-negative real constants a, b, c such that 

(3) ||z|| ^ a + bia{z)y , 
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for every z (z L . Then obviously b,c > , and we say also that a 
is minimal increasing {in any order) not lesser than 1/c. 

Let /i be a measure on / such that intervals are //-measurable sets. 
Functions on / are henceforth read modulo fi (viz. as equivalence classes of 
functions that differ at most on sets of /i- measure zero). Also, let a € M"*^ 
be a delimiting functional on £ ( increasing not lesser than any positive 
order 1/c), and assume that q o / is /i— summable over finite intervals in 
/ for every / € £ . The latter is, for instance, the case if a{z) = || -z || ^ for 
any real p^ 1 and £ C L'^{I, H) which is the metric completion of 



{/ e^' I (/ll/ird^)^ ^ oo} 



i.e. £ is a subset of the Banach space with norm || / ||p : = ( Jj \\ f \\^ dfi) p . 
Definition 4. A causal function F : £ ^ £ is called a-passive, iff 

(4) / a{Ffis))dt,is) ^ f a{f{s))dfiis) , 

Js<t Js<t 

for every / G £ and t € /. 

Remark. If a = ||...||^ for any real p ^ I and || / ||^ is /i-summable over 
/ , i.e. f eLP{I, L), then @ clearly implies Ff G LP{ I, L ) and 

ll-F'f II =^ II f II ■ 

M •> \\p ^- W tf \\p 

Hence, every ||...||P-passive causal function F defines a contraction operator 
on£ n LP{I, L). 

Assume furthermore that r G M_(_ and let on / : = M the measure fi be 
concentrated in {kr \ k (z 7^} with uniform weight fi ( {kr} ) = t , k (z 7^ . 
Alternatively, let / : = {kr \ k ^ 7} with 'the same' measure p, . (Vir- 
tually we deal of course with that discrete situation, even in working on 
the real axis with functions that are constant over intervals [kr, (A; + 1) r) ; 
A; G Z . We retain the integral formalism for simplicity, but the reader may 
optionally re- write the following integrals as sums.) For every / G £ let 
f{t — t) G £, i.e. £ is closed under time shifts by negative integer multiples 
of T. 

For arbitrary A^ G N+ and exciting function e G £ C XL with support on 
[ 0, Nt ) C / , the following holds 

Theorem 1. [Stability of the iterated passive causal process] 

For every a-passive causal function F : £ ^ £ and e (^ £, as stated, 

if g £ L is a process such that g{t) = for t ^ and 

(5) 9it+T) = F[e + g]it) 

for t = nr; n G N, then g is uniquely defined {modulo /i), and for every 
t ^ Nt holds 

(6) \\g{t)\\ ^ a + { - a{e + g)-a{g) dfi )" 

^ J[0,Nt) 

with all constants a,b,c that satisfy ©. Hence g is stable. 
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Remark. A process g ^ ^^ that is recursively generated according to Q by 
iteration of an a-passive causal function F is called an a-passive process. 

The theorem ensures thus that for any excitation of finite duration (and 
with no further restrictions) the a-passive process is necessarily stable and 
in this sense unconditionally stable. 

Note that existence of such a process g is not a priori guaranteed, since 
this obviously depends on the condition that with sq : = e E £ 
also the following functions belong to £ 

(7) s„ := [e + F[s„_i]]<_ G £; for < n < A^ , 

which has to be checked if need be. 

Clearly, condition Q is always true (hence g exists) if £ is a linear space. 
We do not universally premise this, in order to apply the theorem also to non- 
linear situations, where conditions may only be satisfied for sufficiently 
small excitations e . 

Corollary. 

(i) In the special case a = ||...|| estimates holds with a = 0, 6=1, 
c = 1. Then the triangle inequality applies to the integrand of ^ 
and validates the bound 

Wait) II ^ - / \\e\\ dfi . 

■'" J[0,Nt) 

(ii) If N = 1 , i.e. e{t) is a DiRAC excitation concentrated on [0,r) 
(where g = 0), then @ reads simply 

11^(011 ^ a + (- / a{e) dfi )^ , 

^ J[0,r) 

provided that a ( ) = ( which is the normal case ) . 

Proof. Clearly, g ^ L is uniquely defined by the given recurrence relations, 
since e + g { and hence g ) bX the right hand side of © is evaluated only 
up to time t = nr in virtue of the causality of F . 

Furthermore, if A^ ^ n and a ^ || g (nr ) || with any a satisfying Q, then 
with pertinent b, c that satisfy Q 

^ (^(llff(nr)ll -a))i/'= ^ aiginr)) 

= -( / a{g)dii{s) - / a{g)dn{.s)) 

' J s<{n+l)T J s<nT 

= -(/ a{F[e + g])dfi{s) - a{g)dfi{s)) 

' J s<nr J s<nT 

<• recursion formula (JSJ 

< - ( / a{e + g) dfi{s) - a{g) dfi{s) ) 

' J s<nr J s<nT 

<• since F is passive 
= - / a{e + g) - a{g) dfi 

T J[0,Nt) 

T since e{t) = if i^ [0, Nt) . 
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Thus, estimates © applies in the case a ^ ||5'(?ir) || and triviahy other- 
wise. It follows that II g \\ is uniformly bounded on / \ [0, Nt) , hence also 
on / , which is to say that g is stable. D 

5. Stable DSC processes 

In the following, DSC schemes are represented as paired a-passive pro- 
cesses such as dealt with in Theorem ^ and which hence are stable. 

Proposition 3. For every r,s £ J and Tq : /(t) i— > f{t + q) denoting the 
shift operator on L^ ^"^ , 

\ F 

if r -|- s ^ and < J^ is a-passive on 

f TroF„oTs . 

then < ^ J^ ^ IS a-passive on 
\^TrO F^ oTs 

The statement is a direct consequence of Definition |1] D 

Definition 5. 

„, reflection map D? . „ , . .a. \ F„ . 

ihe < . ^ IS called a-passive, iff < ^"- is a-passive 

1 connection map C \ r^ 

in the sense of Definition |^ 

Theorem 2. With every time limited excitation, the DSC process generated 
by a-passive reflection and connection maps is uniformly bounded, hence 
stable. 

Proof. It is sufficient to show that every finitely excited DSC process which 
is generated by a-passive 01 and S can be written as a pair of processes, 
either of which satisfies Theorem ^ 

With H := I U J = {A;r/2 | A; G Z} and the measure on H inherited 
from / and J jointly, the space of all DSC processes is 

£ : = {z={ Zin, Zout ) G {^^) ^ I zUC^k + 1) r/2 ) = z,„( kr ) 
and Zouti kr ) = ZoutiC^k - 1) t/2 ) , for every k eZ} , 

i.e. £ consists of all pairs of functions z = {zin , Zout ) £ (^^)^ > such that 
Zin : H ^f L and Zo^t : H ^> L , respectively, 'switch' at even and odd 
integer multiples of half the time step. So, there is a natural bijection 

£ ^ ^^ X ^■^ 

Z — \ Zin ) ZquI ) ' ^ (, Zin I -' i Zq^i [ 'J ) ) 

in virtue of which the first and second components in £ can be naturally 
identified with L^ and £'' , respectively. 

From section |31 (Definition!^ we recall that for any incident function e € 
L^ X { } C £ supported on a finite interval [ , Nt ) C I , the DSC process 
excited with e and generated by 01 and C is the well-defined function 
2 G (L^) such that z{t) = for t ^ and recursively for < t £ H 

^ \{T_^F^T_^[zout]{t), Zoutit)) ate J. 
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In fact, z is uniquely defined by these relations which provide separate 
recurrence relations for the two processes Zin , Zout \ for instance for z^ and 
t €l 

r^ T ~ T 

Zin{t + t) = Zin{t + -^ ) with t : = t + - € J 

= T.rF^T^ [zout]^^^{t^) by® 



(9) 



By causality Zin enters the last expression only up to argument t , hence (jHl 
is a well defined recursion formula for Zin{t) , t (^ I. Moreover, since prod- 
ucts of a-passive operators are a-passive, F^ T_t/2 Pr ^-t/2 is a-passive 
in virtue of Proposition |21 Hence Theorem ^ applies to Zin which thus is 
stable, just as then also is Zout = F^ T_^/2 [zin + e] ■ D 

6. Model equations 

The typical structure of the DSC algorithm has been outlined in the last 
section. Here we deal shortly with the dynamical model equations - that 
the algorithm has ultimately to solve, and from which the connection and 
reflection maps S, 3? are derived. Most fundamentally, only total fields 
zP, z'^ - not, however, their incident and outgoing components separately, 
enter the DSC model equations. Accordingly we consider only equations 
between total fields, viz. sums of incident and reflected flelds which are of 
the following so called port and node types 

Definition 6. For t e H = lUJ and z{t) = ( Zi„(t) , Zout{t) ) G £ C (£.^ )^, 
(10) ...(t) + .o.aO =: |^„(,) ,^^^j 

In general, viz. prescinding from certain restrictions that can be inferred 
from a principle of near field interaction |He3j . the DSC model equations 
are of the two types 

(11) r'lziJizP] = and :^[4][^"] = , 

with causal functions J'"', 3'^ that respectively switch on / and J . The 
brackets are here short hand for 'back in time running' sequences of states 
(along with their history ) 

[^(±)] := [^]<t(±^) = (^(*(±^) - ^^))n = o.i,2,... , t e H . 

The r/2 time shifts synchronize node and cell boundary switching in (|llj) 
such that the equations hold on all time intervals [t, t + t) , t £ I or t £ J , 
respectively, in harmony with the switching conventions (|10|) for port and 
node quantities. Of course, time shifts by — t/2 would also bring about 
synchronization. However, the resulting equations would not yield explicit 
schemes, in general, or conflict with the causality requirements for 31 and C. 
Inspection of equations (jllj) shows that 9^" affects only the reflection cycle, 
while S'P applies only to the connection cycle. 
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The model equations are treated in more detail in |He3j . In the next 
section the derivation of DSC model equations is exemplified with a heat 
propagation scheme. 

7. A HEAT PROPAGATION SCHEME IN NON-ORTHOGONAL MESH 

The physical interpretation underlying the following application relates a 
smoothly varying (viz. in time and space continuously differentiable, C^-) 
temperature field T, evaluated as T^ at the face centre points and as T" in 
the nodes of a of non-orthogonal hexahedral mesh, to total states z^'" of a 
DSC model. In fact, we derive the model equations for the connection and 
reflection cycles of a DSC heat propagation (diffusion) scheme. Since the 
equations are linear and of dynamical order and 1, respectively - as will be 
seen - they can be processed, following the guidelines of the last section. In 
the end, we display some computational results of a dispersion test carried 
out with this model. 

In order to simplify the notation we follow Einstein's convention to sum 
up over identical right-hand (!) sub and superscripts within all terms where 
such are present (summation is not carried out over any index that also 
appears somewhere at the left-hand side of a pertinent symbol - thus, in 
{—l)'^a'^bx kC the sum is made over A but not over k). 






'U 


/ 




iky 


-». b 




/^f\ 




y^ 


i^ 


7 
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Figure 1. Non-orthogonal hexahedral mesh cell. 
(a) Edge vectors. (b) Node vectors. 

The shape of a hexahedral cell is completely determined by its 12 edge 
vectors {u&)u=o,...,ii- Also, with the labelling scheme of fig^, node vectors 
(^6)^=0,1,2 and face vectors (t/)t=o,...,5 are defined as 



(12) 



and 



J 



1 ^-^3 

4 ^u=o 

i-iy 



(4m+^) 



V (8+20^ "'"(9+2(t+(-l)')) 
'^ \ (4+21,)^ ' (5+21.) ^ ) 



e) A 



^i = 0,l,2 



t = 0,...,5 



(edge vector indices cyclic modulo 12 and A denoting the wedge ('cross') 
product in M^). 

At every face l S {0, ...,5} of a mesh cell and for any given r € M+ 
the following time shifted finite temperature differences in directions ^^6 
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Figure 2. Face vectors and temperature points (nodal section). 



(fi = 0,1,2) form a vector valued function 



(13) ,v^r^(t,r; 



2 (-!)'( T" I 
( 



2fi.+l 



J^P 



t-T/2 



2^ 



.TP 



TP ) 



if/" = ^/2] 
if/u / [^2] 



( [ X ] denotes the integer part of x G M ) . The time increments are chosen to 
attain technical consistence with the updating conventions of DSC schemes. 
They do not destroy convergence, as easily seen: In fact, in the centre 
point of face l the vector iV^T approximates in the first order of the time 
increment r, and of the linear cell extension, the scalar products of the node 
vectors with the temperature gradient VT . Let, precisely, for a fixed centre 
point on face l and e G M+ the e-scaled cell have edge vectors te~ : = e ^e . 
Let also ^V^ T^ denote function ()13() for the e-scaled cell (with node vectors 



6^6 ). Then at the fixed point holds 



(14) 



< ,,6,grad(r) > = ^b-VT 



lim lim — ,V 



B 



^fj. ; 



as immediately follows from the required C^-smoothness of the temperature 
field T. 

To recover, in the same sense and order of approximation, the gradient 
VT from (|1.S|) . observe that for every orthonormal basis {uu)u=o,...,m-i of 
V = M^or C"* , and for an arbitrary basis (/x&)/x=o,...,m-i with coordinate 



matrix (3i, 
are 

(15) 



^ u"^ ; fj, 



b >, the scalar products of any vector a ^ V with nb 



< ijb, a > 



E 



m— 1 
i/=0 



< nb, yU > < yU, a > 

V '^ V ' 



/3>^ 



(At the right-hand side, and henceforth, we follow Einstein's convention to 
sum up over identical sub and superscripts within terms where such are 
present). It follows that 



(16) 



Oiv 






with 



(7.^) ■■=m)T' 
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In other words, the scalar products of any vector with the basis vectors 
^b transform into the coordinates of that vector with respect to an or- 
thonormal basis ,yU via inultiphcation by matrix 7 = {P*)~^ , where Pu = 
< uU , fib > , i.e. f] is the matrix of the coordinate (column) vectors ^6 with 
respect to the given ON-basis ^u , and 7 is the adjoint inverse of f3. 

This applied to the node vector basis ^6 and (|14|) yields the approximate 
temperature gradient at face l as 

(17) yr, = 7^ ,v^r^. 

The heat current into the cell through face t with face vector components 
tf = < J , uu > , u e {0, 1, 2} , is then 

.^gx J = XhJ- NT = XhjTj^ N'^Tf, = ,s^ ,V% , 

Xh denoting the heat conductivity in the cell. 

The heat current through every interface is conserved, i.e. between any 
two adjacent cells C, x with the common face labelled i in cell C, and k in x 
applies 

(19) \J = -IJ. 

Also, the nodal temperature change in cell C is 

(20) ^ = Ajis + y' U), 

where c^ denotes the heat capacity (per volume), V the cell volume, and S 
any heat source(s) in the cell. 

We finally introduce quantities ^zf,,'" {i = 0, ...,5; fj, = 0, 1,2), which still 
smoothly vary in time with the temperature T (and that are hence not yet 
DSC states, but will be later updated as such) 



(21) .z;it) 

and 

(22) .<(t) 



2(-l)'T"|, if^ = W2] 

{2ti+iTP - 2fjiTP)\t^r/2 else 

2{-iy JP\, if/i = [./2] 

,z;j(t-r/2) else. 



From (UHl [THl) follows 

J\t+r/2 = rS^{al\t-2{-ir6^;,'^KTP\t+r/2) 

= .S>^{.z';\t-5TKzl\t+r/2) . 

This, together with (fT9|) . (|22|) . and continuity of the temperature at the 
interface, ^^T^ = ^T^ , implies 



(24) Iz^ 



c 

"' '^l I t+r/2 



'^'^\t + >- 



if^l = [i/2] 



^ U]l\t else, 

which form a complete set of recurrence relations for z^ (given z"). So 
they can be taken as model equations for the connection cycle of a DSC 
algorithm. 
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Equations H20|) discretely integrated in a time balanced form with incre- 
ment T yield 

(25) r''|i+r/2 = '^''\t-T/2 +-^{S + ^^^^t,J\t), 

i.e., with (PT1[^05|) . the recurrence relations 
(26) 

^^ [t/2] I t-T/2 + 2c„y i -^ + 
.^;:|t+./2= I + Ef=0.^'^(^^"|i-r/2-5i'/".^^|t)} if/^=[V2] 

, -^(2^+12^/^ +2/.z^)|t else, 

which provide a complete set of model equations for the reflection cycle 
of a DSC algorithm. Note that the first line, modulo the factor 2(— 1)'', 
always updates the nodal temperature. Of course, this has to be carried out 
only once per cell and iteration cycle. In this - typical - example the dual 
state space concept of DSC ( needed by Johns' cycle ) creates a redundancy, 
which can yet be exploited for process parallelization within either parts of 
the connection-reflection cycle. 

Equations (^H EH]) can be directly taken as updating relations for total 
quantities of a DSC scheme. Alternatively, they may be further processed 
in deriving reflection and connection maps, along with stability bounds for 
the time step. The proceeding is canonical and amounts in essence to a 



straightforward transcription of the model equations along the lines of [Hg 
section 4 (Theorem 1 and corollaries). 

It is quite easy to couple this heat conduction model - within one and the 
same mesh - to a Maxwell field TLM model in the non-orthogonal setting 
|Helj . In fact, with the node vector definition in |Helj . the total node 
voltages are just the scalar products of ^b with the electric field, hence the 
dielectric losses and heat sources per cell are 

(27) S = \aVE^E; = ^^^EJ^^^mI'' 

for a frequency domain (complex) TLM algorithm, cf . |IIe2j : a = 27r/ etan((5) 
denotes the effective loss current conductivity at frequency / in a mesh cell 
of absolute permittivity e and dielectric loss factor tan((5) ; 7 = (/3*)~^ as 
in (|16)) . In Spinner's Maxwell field solver the model couples, in addition, 
to magnetic and skin effect losses. 

Fig El displays the result of a dispersion test, computed in a square mesh 
using non-orthogonal cells. It turns out that the heat conduction properties 
of the mesh are highly insensitive to cell shape and orientation (as of course 
should be the case). 

8. Conclusions 

DSC schemes extend the Transmission Line Matrix numerical method 
in preserving its essential assets. So, they remain stable under quite gen- 
eral circumstances (made tangible with the notion of a-passivity) and are 
amenable to deflection techniques which generalize the stub-loading feature 
of TLM. DSC schemes offer enhanced modeling potentiality, and a wide field 
of future research. 
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Figure 3. Transverse heat conduction over a square mesh 
using non-orthogonal cells. A Heaviside temperature step is 
imposed on one side and the transient temperature computed 
at the opposite side, assuming adiabatic boundary conditions 
on all but the heated sides. 

DSC results (+) are plotted over analytical solution (curve). 
(a) The mesh, (b) Horizontal (c) vertical propagation. 
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